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Abstract 

We connect two different extensions of Boltzmann's kinetic theory by requiring the same station- 
ary solution. Non-extensive statistics can be produced by either using corresponding collision rates 
nonlinear in the one-particle densities or equivalently by using nontrivial energy composition rules 
in the energy conservation constraint part. Direct transformation formulas between key functions 
of the two approaches are given. 
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Power-law distributions in nature are nearly as common as Gaussian distributions, which 
is the limiting distribution emerging for the properly scaled sum of infinitely many inde- 
pendent random variables. The power-law tailed distributions have escaped somehow the 
strength of the central limit theorem: either by being composed from long-tailed individ- 
ual distributions, or by featuring finite-size effects in a general sense. To the latter cat- 
egory belong unscreened long-range forces and long-time memory effects; to the former a 
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(mult i) fractal phase space occupation. Since power-law tailed distributions can 
many areas from particle physics and astrophysics to financial market models 
righteously suspected that these belong to some more or less universal stationary states of 
complex dynamicses. The question arises whether there are characteristic common features 
in these dynamicses and what they then are. 



Non-extensive thermodynamics has been developec 



in the past two decades statistical 



theory to deal with such physical stationary statesjj, ^,|f3]. Initially based on mathematical 
investigations of a generalized definition of Boltzmann's entropy, a never-decreasing macro- 
scopical state-parameter intimately connected to statistical probabilities of microstates, this 



theory soon started to study dynamicses possibly 



eading to such states. Studies of anoma- 



B,y,3 



revealed that either a nonlinearity 



lous diffusion, random walk, and noisy equations [ 
in the one-particle probability or a nontrivial interaction with the environment, which de- 
pends on the observed low degree of freedom subsystem itself, may be responsible for such 
a nontraditional behavior. 

While the above approaches all mound in the study of an extended Fokker-Planck prob- 
lem, the set of all dynamical evolutions leading to a thermodynamical state is much wider. 
In general a considerable fluctuation of the intensive parameters may establish a non- 
Boltzmannian canonical distribution Q]. Another classical field of establishing thermody- 
namics (and hence also non-extensive thermodynamics) is kinetic theory in general. The 
basic assumption of Boltzmann, that the time evolution of many systems has a fast, micro- 
reversible and ergodizing component, which he comprised into a collision integral, opens up 
studies of more complex non-linear dynamical theories. 

Boltzmann's original kinetic theory is based upon a collision rate, which is multilinear in 
the one-particle densities and symmetric both in the colliding partners and against time re- 
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versal. Based on these properties the H-theorem was derived, funding the quantity entropy, 
and giving a way towards a microscopic establishment of thermodynamics. The Boltzmann 
entropy is an extensive (additive) quantity reflecting the independence of microstate proba- 
bilities in weakly interacting subsystems. The classical Boltzmann equation restricts possible 
stationary distributions to satisfy a product rule, once an addition rule for energy is given. 

Recent generalizations of the Boltzmann equation hook in exactly at this point: either 
the product formula is generalized^,^!, interpretable as mimicking nontrivial many-body 
correlations in the equilibrium state, or the addition of energy is replaced by a more general 
formula accounting for an in- medium interaction energy shared by the colliding part- 
ners. These two approaches may lead to the same, non-standard stationary one-particle 
distribution. The aim of this work is to clarify the interrelations between them and to 
present corresponding mathematical formulas. 



In Ref . |H| possible extensions with nonlinear dependence of the collision rate on the one- 
particle densities have been considered, this approach we shall refer here as the nonlinear 
Boltzmann equation (NLBE). The rate in this approach consists of production and blocking 
factors, both in gain and loss terms. Considering a simple 1 + 2 3 + 4 two to two body 
collision, the rate of change of the one-particle phase space density, fi = f(pi), is given by 

h = ^1234 (a 3 bxa 4 b 2 - a^a^) (1) 

J 234 

with cij = a(fi) being the general production and 6j = b(fi) the blocking factors. The 
transition probability rate factor, W1234 may contain a 1—3 and a 2—4 symmetric contribution 
also from correlations between production and blocking. This is a general extension of the 
Boltzmann and the Boltzmann Uhling-Uehlenbeck equations with respect to nonlinearity 
of the collision rate. The standard theories are recovered for a(f) = f and b(f) = 1, or 
b(f) = 1 ± / respectively. In this case the stationary distribution is governed by the ratio 
n(f) = a(f)/b(f) which becomes the traditional Boltzmann factor: 

K(f eq ) = eM-E/T). (2) 

This result assumes that in two-body collisions momenta and energy are composed additively 
(E 1 + E 2 = E 3 + £4, pi + p 2 = P3 + f>i). An H-theorem can be proven with the generalized 
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expression for the entropy, 

Sk = I <r(/i) (3) 

where the integration is over the one-particle phase space. It turns out that cr(f) is related 
to the previous quantities via 

a'{f) = - In «(/). (4) 

Recently another approach to generalize Boltzmann's original treatment has been 
proposed^]. Here the (multi) linearity of the collision rate is kept, blocking factors are 
not applied, but the additivity of the energy during the micro-collisions is replaced by a 
more general requirement: only a given function of the individual energies, physically stand- 
ing for the total two-particle energy, is conserved, h(E\,E2) = h(E$,E/±). The function 
h(x, y) describes a general, non-extensive energy composition rule for the two-body system. 
If this is chosen with the property of associativity, h(h(x, y),z) = h(x, h(y, z)), then its most 
general form is related to a strict monotonic function, X(x): 

h(x,y) = X- 1 (X(x) + X(y)) . (5) 

The function X(x) is a mapping of the non-extensive composition rule to the addition rule, 
it is unique up to a real multiplicative factor. The stationary distribution in this case is 
given by 

f eq = exp(-X(E)/T). (6) 

The H-theorem can be proven for Boltzmann's original construction, but this quantity may 
be interpreted as the total of the additive mappings of non-extensive entropy contributions: 

S B = X S (S N ) = -//In/. (7) 

By requiring a connection to the stationary state described by the previous approach, it 
turns out that the same composition rule should be applied to the energy (scaled by the 
temperature) and to the entropy, i.e. X(E)/T = X S (E/T) with the respective mapping 
function for energy and entropy. The non-extensive Boltzmann equation (NEBE) approach 
is still multilinear in the one-particle distributions. In parton cascade simulations, by ap- 
plying the above rules, the non-extensively composed total energy, E to t = X~ x i^ i X{Ei)), 
is conserved. 



4 



Of course, one may consider to apply both extensions, the nonlinear density dependence 
and the non-extensive energy composition rule. In order to support a given non-extensive 
thermodynamics, however, each one alone suffices. In this respect these two generalizations 
are equivalent and transformation formulas can be obtained between them. In the followings 
we review this correspondence. 

A common stationary distribution relates the n(f) and the X S (E/T) functions: 

f eq = n-\e- x ) = e- x 'M (8) 

with the argument x = E/T. This allows for obtaining both n(f) or X S (E/T) to a given 
distribution, or converting these two basic functions into each other. Based on this the 
non-extensive entropy formulas © and (J7J) can be derived, the relation 

a'(f) = -\nK(f) = X; 1 (-\nf) (9) 

is just a consequence thereof. Applying this relation to the Boltzmann entropy (JJJ) we arrive 
at 

X S (S N ) = I fX s (s(f)), (10) 



with the one-particle contribution s(f) = cr'(f) at the generalization of the mapping of the 
non-extensive entropy to an additive entropy measure, Sb- The basic entropy transformers 
of the two approaches, <r(f) and X S (E/T) must be related by © in order to describe the 
same non-extensive thermodynamics in the stationary state. As a trivial consequence the 
same composition rule, h(x,y), is applied to the energy and to the entropy up to a scale 
factor T. 

It is noteworthy that the Tsallis entropy can actually be obtained by considering St = 
I f s (f)- This very expression fulfills an H-theorem in the case of Tsallis distribution, the 
corresponding Tsallis and NLBE a'(f)-s being linear functions of each other (see below). 
In the general case they are, however, different. The deformed exponential and logarithm 
functions, considered in several approaches to the non-extensive thermodynamics, are simply 
related to the scaled mapping function l^J : 

ex Pdef (t) = e -*-(-*>, \n def (f) = -X;\-lnf). (11) 
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Finally we demonstrate the above relations for some well known stationary distributions, 
relating them to specific non-extensive composition rules, production and blocking factors, 
or expressions for the total entropy, respectively. 



The Maxwell-Boltzmann distribution, or in general the Boltzmann-Gibbs distribution, 
f(E) = exp(—E/T) is recovered as a particular case in both approaches. There is no 
blocking factor, b(f) = 1, the production factor is linear, a(f) = f. The mapping function 
is the identity, X s (t) = t and the energy (and entropy) composition rule is the simple 
addition, h(x, y) = x + y. 

The Tsallis distribution, f(E) = (1 + a£)" 1/aT arises for X s (x) = - ln(l + arc) [15[ by the 
energy composition rule h(x, y) = x + y + axy (|. The production to blocking factor ratio 
becomes n(f) = a(f)/b(f) = exp((l — f~ aT )/aT). The Tsallis index is given by q = 1 — aT. 
Note that if the parameter a in the generalized composition rule is temperature indepen- 
dent, then the resulting Tsallis index will depend on temperature. The NLBE entropy is 
given by S K = jz^Kf/q ~ f), the Tsallis entropy by S T = j^Kf 1 ~ /)■ This is not a 
principal difference, because whenever an H-theorem is fulfilled with a given cr'(f), a linear 
combination, Aa'(f) + B also fulfills the same H-theorem. The choice A = q and B = 1 
transforms Sk into St in this case. 

Its mapping to an additive quantity coincides with the Renyi entropy jig], Sr = X s (St) = 
In / f q , when the distribution / is normalized to one. 

A variant of the Levy distribution, f(E) = exp(— (E/T) v ), with a fractional power v 
can be achieved by using the production to blocking ratio = exp(— (— In f) 1 ^) in the 
NLBE. The mapping to additive quantity is done by X{x) = x v , generating the abstract 
composition rule: h(x,y) = (x v + y v ) l ^ v '. 

A pure power-law distribution, f(E) = (bE /T)~ , may be produced by using «;(/) = 
exp(— f q ~ x /b) or equivalently by mapping the energy to an additive quantity via X(x) = 
Inbx. The corresponding composition rule is given by h(x, y) = bxy. The NLBE entropy 
becomes Sk = I f q /bq, the Tsallis entropy is given by St = J f q /b in this case. For 6 = 1 
they coincide. 

The so called quon distribution, f(E) = eE /r + (including the Fermi (q — 1) and Bose 
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(q = — 1) distributions as particular cases), is achieved by the production to blocking ratio: 
K (f) = jz^j (the choice a(f) = f and b(f) = 1 — qf leads to the Uhling-Uehlenbeck 
equation). The corresponding mapping to additive quantity is given by X(x) = hx(q + e x ) — 
ln(g + 1) (X(0) = is achieved by a subtraction of a finite constant. This does not work for 
q — — 1 signaling the divergence due to Bose condensation.). The generalized addition rule is 
somewhat complicated: h(x, y) = x + y + In (l + q(e~ x + e~ y ) + q(q — l)e~^ x+y ^ — ln(g + 1) 
features a "Pauli potential", a pair energy leading to an exact Fermi distribution. In the 
NEBE approach this can be achieved without considering a third (blocking) particle at each 
collision. The density for the Sk entropy is given by a(f) = —f In / — (1/q — f) ln(l — qf). 

Another distribution, proposed in Ref. 3] for applying to a kinetic theory of relativistic 
particles, applies the mapping function X(x) = \ ■ asinh(fcc). Another long-known formula, 
the mapping of relativistic velocities to an additive rapidity variable, follows from Einstein 
velocity-composition formula: X(y) = c- atanh(t>/c). 

Some different distributions may be united in two- (or more) parameter classes. For 
example, X(x) = ^ ((1 + ax) q — 1) for small q but arbitrary large x approaches Abe's 
logarithmic formula leading to the Tsallis distribution, while for large a comes close to 
X(x) = x v leading to the Levy distribution jisl . 

Finally non-associative composition rules can be simulated in computerized parton cas- 
cade simulations, too. In this case, however, no mapping can be found to an additive, 
statistical quasi-energy. The individual energies after a micro-collision, although random in 
a certain kinematical range, are no more one-variable functions of the respective energies 
of the incoming pair, but depend on both energies. Non- associative composition rules, on 
the other hand, would not be able to converge in the thermodynamical limit of repeated 
compositions of compositions. Therefore they can probably be only of pure mathematical 
interest . 

We have carried out numerical simulations so far for the Tsallis case^j]. The nonextensive 
Boltzmann equation is simulated by = 10 5 — 10 6 test particles having initially random 
momenta pi in a certain range and evolving by pairwise changes according to the rules 

Pi + P2 = P3 + Pa, X(E( Pl )) + X{E(p 2 )) = X(E(p 3 )) + X(E(p 4 )) (12) 
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massless Tsallis kinematics (E=4, P=3.5) 
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FIG. 1: Kinetic constraint for momentum changes in a pairwise elastic micro-collision using Tsallis 
rules for the energy composition and massless free dispersion relation. The inner curve corresponds 
to a = 0.125, the middle one to a = and the outer one to a = —0.125. 
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FIG. 2: Tsallis distributions obtained numerically by solving the corresponding nonextensive 
Boltzmann equation. 

with the energy mapping X(E) = -ln(l + aE). This simulation conserves the quasi energy 
X(Etat) = J2iX(E(pi)). Fig|T]shows pairs of momentum vectors satisfying the above condi- 
tions for massless (extremely relativistic) particles, i.e. ior E(pi) = \pi\. The two-dimensional 
section of the surface of momenta is an ellipse for a = (the traditional Boltzmann case) 
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only, for positive or negative values of this parameter a fourth order curve is drawn. Final, 
stationary distributions of the one-particle bare energy, f(E) with E = \p\ are then Tsallis 
distributions. FigEI shows examples of these distributions. These are numerical demonstra- 
tions of how to achieve a Tsallis distribution from arbitrary initial distributions. 

In conclusion we have demonstrated that quite different ansatze to generalize Boltzmann's 
classical kinetic theory produce the same stationary distributions of non-extensive thermo- 
dynamics. Upon this equivalence there always exist a strict monotonous mapping to an 
additive entropy from more exotic entropy definitions, as long as the energy composition 
rule applied in the micro-events (collisions) are associative. It is natural to assume that the 
relevant rules in the thermodynamical limit have to become associative. The H-theorem 
holds. Several known non-Boltzmann distributions belong to a physically intriguing non- 
extensive energy addition rule. This gives hope to find realizations in nature by studying 
the pair interaction mechanism and its thermodynamical limit. Nonlinear production and 
blocking factors may formally be replaced by non-extensive energy formulas if only the sta- 
tionary distribution is asked for. The second method, not using phase space blocking factors, 
allows for computer simulations with a resource effort not worse than 0(N 2 ). 
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